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Abstract 



! ! Recently, the concept of rotational mixing has been challenged by some authors 

(e.g. Hunter et al. 2008). We show that the excess N/H is a multivariate func- 
Qh' tion f{M, age, vsvai, multiplicity, Z). To find a correlation of a multivariate 
q , function with some parameter, it is evidently necessary to limit as much as pos- 
' sible the range of the other involved parameters. When this is done, the concept 
, of rotational mixing is supported by the observations. We also show that the 
sample data are not free from several biases. A fraction of ~ 20 % of the stars 
may escape to the relation as a result of binary evolution. 



f^. 1 Introduction 

m ! 

• In the late 70', the reality of mass loss in massive stars was debated as well as its 
\ effect on the evolution. When this became accepted at last (Chiosi & Maeder 
C^) • 1986), the next question was whether the large differences in the populations 
of massive stars (for example the WR/O and WN/WC star numbers) were due 
to differences in mass loss with metallicity Z. Today, the debate concerns the 
reality of rotational mixing and its effects on massive star evolution. History 
never reproduces itself similarly, nevertheless there is a great parallelism in these 
debates, which are normal steps in the progress of knowledge and finally lead 
jJJ ■ to a better understanding of stellar physics and evolution. 

There is an impressive list of consequences of stellar rotation fMaeder &i 
Meynet 2000), many of which are supported by observations: about the stellar 
shape, the temperature distribution at the surface, the mass loss and its asym- 
metries, on the size of the cores, the tracks in the HR diagram, the lifetimes, 
the surface composition, the chemical yields, the ratios of the different kinds of 
massive stars (blue, red supergiants, WR stars), the types of supernovae and 
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Figure 1: The differences in log(N/H) as a function of the initial masses at 3 stages 
during the MS phase for models with Z = 0.02 with average rotation velocities of 
217, 209, 197, 183, 172, 168 km s _1 for respectively 12, 15, 20, 25 40 and 60 M . 
These 3 stages are indicated by the value of the central H-content X c . 

the remnant masses, etc. In this Liege Colloquium, the reality of the rotational 
mixing was disputed by several authors on the basis of new observations, in 
particular the VLT-Flames survey (Hunter et al. 2008). Thus, we concentrate 
here on this problem, firstly by recalling some theoretical predictions concerning 
surface enrichments and secondly by carefully examining the observations. 

Recall of theoretical predictions concerning rotational mixing 



In the mass range of ^10 to 20 M Q considered here, mass loss has a 
limited importance during the Main Sequence (MS) phase. The changes of 
abundances are expected to be mainly due to rotational mixing. The main 
effect producing element mixing is the diffusion by shear turbulence, which 
itself results from the internal SI gradients built during evolution. To a smaller 
extent, meridional circulation makes some transport, however mainly of angular 
momentum. Mixing brings to the surface the products of CNO burning: mainly 
14 N and 13 C enrichments, 12 C is depleted with limited 4 He enrichment and 16 
depletion. 
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Figure 2: The evolution of the differences in log(N/H) duringthe MS phase as function 
of the actual rotation velocities (sini being equal to 1 here) for models of 20 M© 
with Z — 0.02 and different initial velocities. 

Fig. 1 shows the predicted variations of \og(N/H) during MS evolution as 
a function of the initial masses (Meynet & Maeder 2000). (N/H) is here the 
abundance ratio of N and H in numbers (the relative differences in mass and 
number are the same). Without rotational mixing, there would be no enrich- 
ment until the red supergiant stage. Rotation produces an increase (depending 
on velocity v) of N/H during the MS phase. The N excesses also depend on the 
ages t. The increase is modest during the first third of the MS phase, because 
the elements need some time to reach the surface, then it is more rapid. The 
N enrichments are larger for larger masses M. Thus we see that the N excesses 
are multivariate functions 

A\og(N/H) = f(M, t, vsini, multiplicity, Z) . (1) 

Models with lower initial metallicities Z have higher N enrichments for given 
M and v (Maeder & Meynet 2001; Meynet et al. 2006). The excesses become 
very strong at metallicities Z as low as 10~ 8 . 

Fig. 2 shows, for given M and Z , the evolution of the N excesses with age 
for different initial velocities. This figure shows that there can be no single 
relation between vsini and the N/H excesses for a mixture of ages. The same 
is true for a mixture of different masses. 

Binarity may also affect the N and He enrichments due to tidal mixing and 
mass transfer. A binary star with low rotation may have a high N/H due to 
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tidal mixing or due to the transfer of the enriched envelope of a red giant. At 
the opposite, a binary star may also have a high vsini and no N/H excess, 
in the case of the accretion of an unevolved envelope bringing a lot of angular 
momentum. A nice illustration has been given (Martin, this meeting). 

Study of the observed N/H excesses 
The M, age and Z dependences 

Let us start by examining the mass, age and Z dependences of the N/H excesses. 
The data for different groups of stars with different Z are summarized in Table 
1. In the Galaxy (Z w 0.02), the main recent data sources (Herrero 2003; Venn 
& Przybilla 2003; Lyubimkov et al. 2004; Huang & Gies 2006; Trundle et al. 
2007) support significant excesses of He or of (N/H). In the lowest mass range 
studied (6.6-8.2 M Q ), small excesses of He/H are still present (Lyubimkov et 
al. 2004). In the LMC (Z « 0.008), the excesses are larger (Hunter et al. 2007; 
Trundle et al. 2007). In the SMC (Z « 0.004), still much larger N excesses 
are observed (Venn, Przybilla 2003; Heap & Lanz 2006; Trundle et al. 2007; 
Hunter et al. 2007). 

Table 1: The largest Alog(iV/ H) values observed for different types of stars in the 
Galaxy, LMC and SMC (differences in dex with respect to the local values in the 
considered galaxy). The average is equal to about the half of the indicated values. 

Types of stars Alog(N/H) Galaxy A log(N/H) LMC A log(N/H ) SMC 



O stars 0.8 - 1.0 1.5 - 1.7 

B-dwarfs M < 20 M 0.5 0.7 - 0.9 1.1 

B giants, supg. M < 20 M 1.1 - 1.2 1.5 

B giants, supg. M > 20 M© 0.5 - 0.7 1.3 1.9 



These data show the following facts, consistent with theoretical predictions: 

• On the average, the N enrichments are larger for larger masses. 

• The N enrichments are larger at lower Z. 

• The He and N enrichments increase with the distance to the ZAMS 
(Huang & Geiss 2006). They are even larger in the giant and supergiant 
stages (Venn & Przybilla 2003). This property is also well observable for 
example in Nil (Fig. 34 of Hunter et al. 2007) and in NGC 2004 (Fig. 2 
by Trundle et al. 2007). 
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The vs'mi dependence of the N/H excesses 

Several correlations of the N or He excesses with the observed vs'mi have 
been performed. Huang and Gies (2006) and Lyubimkov et al. (2004) find a 
correlation of the He excesses with vsini for B stars in the upper part of the 
MS band in agreement with model predictions. 

In other comparisons (Hunter et al. 2008; Langer, this meeting; Brott, this 
meeting), the authors conclude that "the observation . . . challenges the concept 
of rotational mixing". They claim that "two groups of core hydrogen burning 
stars . . .stand out as being in conflict with the evolutionary models". Group I 
contains rapid rotators whith little chemical mixing, while Group II consists of 
low rotators with large N enrichments. We clearly disagree with the conclusions 
of Hunter et al. (2008), which mainly result from the fact that, instead of Eqn. 
(1), their analysis implicitly assumes that 

AIog(JV/ff) = /(vsini) . (2) 

We note the following points: 

• Their sample contains a mixture stars in the mass interval of 10 to 30 
M . Fig. 1 shows that over this mass interval, the N/H excesses vary as 
much as by a factor of two for a given rotation velocity. 

• The sample by Hunter et al. consists of stars in extended regions around 
the LMC clusters Nil and NGC 2004. As stated by the authors, their 
sample also contains field stars, which do not necessarily have the same 
age or degree of evolution as the cluster stars. Thus, large differences 
of N/H are possible for given M and vs'mi. Also, the two clusters do 
not have the same ages, Nil being younger than NGC 2004, so that the 
stars near the turnoff of Nil have a mass of about 20 M Q , while this is 
about 14 M Q for NGC 2004 according to the HR diagrams by respectively 
Hunter et al. (2007) and Trundle et al. (2007). 

• The completeness of the binary search is unknown. 

To limit the severe effects of mass and age differences, we consider separately 
the two clusters. In Nil, we limit the sample to the stars in the mass range 14 
to 20 M Q on the basis of the data provided by Hunter (2008) and in the formal 
MS band as given by Fig. 34 from Hunter et al. (2007). In NGC 2004, we take 
the mass interval 13 to 16 M Q (same source) and in the formal MS band from 
Fig. 2 by Trundle et al. (2007). Figs. 3 and 4 show the results. Ideally these 
mass intervals should even be smaller. 

For Nil, we see that Group I (stars with high vsini and low N/H) has 
essentially disappeared. There remain only 2 evolved binary stars, which is 
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Figure 3: The N abundance (in a scale where logff = 12.0) as a function of «sini 
for the MS stars (black dots) in Nil with masses between 14 and 20 Mq according 
to Hunter (2008). The binaries are shown by a square. The evolved stars in a band 
of 0.1 dex in log Teg beyond the end of the MS are shown with open symbols. The 
gray band indicates uncertainties of ±0.25 dex. 

consistent with some scenario of binary evolution. We suspect that Group I 
was also largely formed by stars of smaller ages and/or lower mass stars, where 
for a given vs'mi the N excesses are smaller. A support to the latter possibility 
comes from the fact that the average mass of the stars in the region of Group 
I is 12.8 M (for vsini > 180 km s" 1 and 12 + log(N/H) < 7.30) and 17.1 
M for stars with vsini > 180 km s _1 and 12 + log(N/H) > 7.30. Group II is 
mainly made from evolved stars (open circles), which explains the low velocities 
and high N abundances. The two remaining stars in this group can easily be 
stars with a small sinz, especially more than the sample data by Hunter et al. 
(2008) is largely biased toward low rotators. The ratio of star numbers with 
vs'mi > 250 km s _1 to those with vsini < 100 km s _1 is 0.14, while in the 
clusters studied by Huang & Gies (2006), this ratio amounts to 0.40! Thus, the 
sample by Hunter et al. contains a large excess of slow rotators. We conclude 
that the bulk of stars in Nil shows a relation of the excess of N/H depending 
on vs'mi (the mean square root of the data for the MS band stars is 0.23 
dex from the data by Hunter (2008), the scatter in vsini is not given). The 
amplitude of the (N/H) is about 0.6 dex for velocities of 200 km s _1 , slightly 
higher than at the corresponding mass in Fig. 1 for Z = 0.02. 

For NGC 2004, the results are essentially similar. For most stars there is a 
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Figure 4: The N abundance as a function of vsini for the MS stars in NGC 2004 
with masses between 13 and 16 Mq. Same remarks as for Fig. 3. 



relation between the excess of N/H and vsini. For Group I, there remains only 
star nb. 100 (which is not much for a group). This star is interesting. Its mass 
is 13 M©, vsini is 323 km s _1 , the highest of the whole sample. In reality, the 
velocity is still higher because the authors do not account for gravity darkening 
(Fig. 4 of Hunter al. (2008) tends to support this remark). This star might be 
a re-accelerated binary or simply a younger star in the field. Another possibility 
(which we favour) is that its parameters have been incorrectly appreciated due 
to the extreme rotation. In this respect, a logg vs. logT c fr diagram, e.g. 
Fig. 16 by Meynet & Maeder (2000), shows that a too high mass is assigned 
to a fast rotating star if in the log g vs. logT e g- diagram its mass is determined 
from non rotating models. Thus, the mass of star nb. 100 could be lower than 
13 Mq which is the lower bound of our sample. To know whether this is what 
occurred for this star, the whole reduction process should be redone. For Group 
II, a large fraction consists as for Nil of evolved stars. Again for NGC 2004, 
the N enrichments increase with vs'mi in agreement with theory. 

We also note that the data used in the analysis by Hunter et al. (2008) 
are subject to several biases. Firstly, the sample contains no Be stars, while 
their number fraction is about 15 to 20 % in the LMC. This contributes to bias 
the sample toward low velocities. A second source of bias is that the vsini 
determinations are based on models assuming that the stars are uniformly bright, 
with no account given to gravity darkening (now an observed effect). Thirdly, 
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in the values of gravity used to estimate the masses, no account is given to the 
gravity change due to rotation. The effects of evolution and rotation should be 
disentangled before any mass is assigned. 

Conclusions 

We conclude quite logically that to find a correlation for a multivariate function 
like N/H with some parameter like vsin-i, it is necessary to limit as much as 
possible the range of the other involved parameters. Otherwise, the conclusions 
may be erroneous. We note that data samples limited in mass and ages support 
a N enrichment depending on rotational velocities. Stars beyond the end of the 
MS phase do not obey to such a relation, because their velocities converge 
toward low values (see Fig. 12 by Meynet and Maeder 2000). A fraction, which 
we estimate to be ~ 20 % of the stars, may escape from the relation as a result 
of binary evolution, either by tidal mixing or mass transfer. 
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